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1. Introduction 

Heavy quarks are emerging as valuable probes for the study of quark gluon plasma 
produced in relativistic heavy ion collision. This has its origin in the large mass of 
heavy quarks which lends them quite a few advantages, viz, they are produced at a time 
2^7- which is smaller than the formation time of quark gluon plasma. Their large 
mass ensures that their production can be calculated reliably using perturbative QCD 
and they may not be deflected substantially from their path due to collisions with quarks 
and gluons and due to radiation of gluons. The drag experienced by the heavy quarks 
due to these collisions and radiations however leads to medium modification of their 
production which is quite similar to those for light quarks, as leading particles [TJ [2] . 
Recent calculations which treat the so called 'dead cone' with more care [3], also show 
that heavy quarks lose energy in a manner quite similar to light quarks [I]. 

Recently it has been pointed out that correlations of heavy quarks (charm and 
anti-charm) can add a richness to these studies by adding several features [5] . Consider 
for example heavy quark production at lowest order of pQCD. They would be produced 
back to back. The two members of the correlated pair may in general lose different 
amounts of energy as they may cover differing path lengths in the plasma. However if 
they do not change their directions, they would continue to remain back-to-back. Now 
consider that there is a strong flow and that the heavy quarks take part in flow [6]. It 
is now possible that one of them is produced with a transverse momentum parallel to 
flow velocity and its momentum will increase, while the momentum of the other will 
decrease. In fact if the radial flow velocity vt > Pt/ E (charm), the charm will change its 
direction and the back-to-back correlation may change to forward correlations. When, 
however, the flow velocity is not collinear with the momentum, the final momenta will 
be separated by < < it. Thus while the energy loss is not likely to alter the angular 
correlation of heavy quarks at lowest order in pQCD, a strong elliptic flow will bring 
in some interesting and rich structure, the analysis of which could throw some light on 
interplay of energy loss and flow. 

There is, however, a substantial production of heavy quarks at next to leading order 
in pQCD. The NLO process gg — > QQ can proceed in two ways(among others). Either 
one of the final state gluons in the process gg — > gg splits (g — > QQ) or one of the heavy 
quarks radiates a gluon following gg — > QQ. The pair is expected to be collinear in 
the first case and deviate back-to-back in the second case. The processes where gluon 
is emitted from the external legs will fill up the region < < 7r. Now energy loss 
will alter the correlations in a complex manner. If our assumption on heavy quarks 
not changing direction due to energy loss largely holds then px integrated correlation 
is likely to remain unchanged. However if we now study the correlation for different 
cuts on pt, some interesting patterns may emerge. Different heavy quarks lose different 
momenta! 

We can now discuss correlated decay of charm-anti-charm into electron-positron 
pair. The invariant pair mass distribution of electron pair obtained from decay shows 
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interesting features. It is seen earlier that large suppression of heavy quark as seen 
through Raa-i results in increase of D mesons as well as single electron spectrum at 
low momentum by a few percent. This characteristic increase is quite different from 
enhancement due to Cronin effect [7] and is found to be due to large effective drag 
upon charm by thermalized medium. The invariant pair mass distribution of electron 
pair obtained from decay shows similar feature from effects due to energy loss by 
charm quarks [8]. The electron pairs pile up at low invariant mass region resulting 
in characteristic enhancement in electron distribution. 

In the following we study some of these features of the correlation of heavy quarks 
with collision of lead nuclei at 2760 GeV/nucleon as an example. The paper is organized 
as follows. Sec. 2 contains formalism for charm production cross-section from pp 
collisions and lead on lead collision at LHC energy. Sec. 3 contains an empirical model 
of energy loss employed to determine the medium effect on charm correlation. Sec. 4 
presents our results and discussions on azimuthal correlation and correlated charm decay. 
Finally Sec. 5 gives the summary followed by acknowledgement and bibliography. 



2. Formulation 



The correlation of heavy quarks produced in pp collisions is defined as 

da da 



Ei E 2 



C 



(1) 



d?p x d 3 p 2 dy 1 dy 2 d 2 p T i d?p T2 
where yi and y 2 are the rapidities of heavy quark and anti-quark and pxi are the 
respective momenta. 

At the leading order, the differential cross-section for the charm correlation for 
proton on proton collision is given by 

da 



C 



LO 



dyi dy 2 d 2 p T 
One can now calculate [5J j9] 



<Hpti + PT2) 



(2) 
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where px and y\^ 2 are the momenta and rapidities of produced charm and anti-charm and 
x a and Xb are the fractions of the momenta carried by the partons from their interacting 
parent hadrons. These are given by 



X n 



M T 

7^ 



e yi + e 



!J2 ' 



M T 

7^' 



-m 



+ e 



(4) 



where Mr (= yrriQ +pj<), is the transverse mass of the produced heavy quark. The 
subscripts i and j denote the interacting partons, and f/j are the partonic distribution 
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functions for the nucleons. The invariant amplitude, \M\ 2 in differential cross-section 
da/di is taken from ref. [ID] . 

The processes included for LO calculations are: 

9 + 9 -> c + c 

g + g ->■ c + c . (5) 
At Next-to-Leading order the subprocesses included are as follows: 

9 + 9 c + c + g 
q + q -» c + c + g 

g + q(q) -» c + c + g(g) . (6) 

The eq. [I] gives the correlation of heavy quarks from initial fusion in proton-proton 
collision. The azimuthal correlation of heavy quark for Pb+Pb collision at given impact 
parameter is given by 

E c E 5 dNAA = T AA E C E- c d(Tpp ( 7) 

c d 3 p c d 3 p 5 AA c c d 3 p c d 3 p 5 

For lead on lead collisions at LHC, we have used T A a — 292 fm~ 2 for b = fm 
We have used CTEQ5M structure function. The factorization, renormalization, and 



fragmentation scales are chosen as 2ym 2 + p\ and the charm quark mass, m c has been 
taken as 1.5 GeV. 

3. Energy Loss Formalism 

We use the empirical model for the energy loss for charm quarks proposed in one of our 
earlier paper. We perform a Monte Carlo implementation of our model calculations and 
estimate the azimuthal correlation as well as correlated decay of charm pair with charm 
cross-section determined using NLO-pQCD calculations. 

We assume that the energy loss of heavy quarks proceeds through the momentum 
loss per collision is given by, [TT] 

(4p)< = a (p*f , (8) 
so that one can write 

dV = _^P (g) 
dx A 

where a and (3 are parameters with best values at y/s= 2760 GeV/nucleon taken from 
publication by Younus et al [1] and A is the mean free path of the charm quark, taken 
as 1 fm throughout. Thus the momentum of the charm quark after n collisions will be 
given by 

p n+ l =p n - (Ap) n (10) 
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The probability for the charm quark to have n collisions, while covering the path length 
L is given by 

P(n,L) = ^e^/\ (11) 

So now we estimate the largest number of collisions- N, which the charm quark 
having momentum pt can undergo. Next we sample the number of collisions n, which 
the charm undergoes from the distribution 

p(n) = P(n,L)/J2P(n,L) (12) 

n=l 

to get the final momentum of the charm(anti-charm) quark. 

Next we fragment the charm quark using Peterson fragmentation function given by 
D. We have assumed that D(z), where z = Pd/Pc, is identical for all the D-mesons, [T2] 
and 

z[l - 1/z - e p /(l - z)} 2 
where e p is the Peterson parameter and 

dzD(z) = l. (14) 

We have kept it fixed at e p =0.13. 

Then we have included semileptonic decay of D(D) mesons by parameterizing 
electron distribution function taken from Ref. [13]. Finally we show our results for 

dN C c/dA(j), dN DD /dA(f), E c EcdN/d 3 p c d 3 Pc and dN/dM e+e -. 
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Figure 1. (colour on-line)Comparison of D mesons azimuthal spectrum for two 
different structure functions. 
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Figure 2. (Colour on-line)(i7V/<iA</i vs A(j> of cc pair for (upper left)pr < 2.0 GeV, 
(upper right)p T > 2.0 GeV. (lower left)p T < 6.0 GeV, (lower right )p T > 6.0 GeV. 



4. Results and Discussions 

We have used NLO-MNR code [H] with CTEQ5M structure function for estimating 
charm cross-section for all leading and next-to-leading pQCD processes. The scaling 
factor used is 2\jm 2 c + p\ with m c =1.5 GeV. In this paper we have used two different 
values for parameter '/3'=1.0 for B-H type and /3=0.5 for LPM type of energy loss 
mechanisms respectively Correspondingly a = 0.12 for B-H type and 0.25 GeV 1 / 2 
for LPM type are taken as the best values at y/s = 2760 GeV/nucleon. The entire 
calculation is done for central collision (b=0fm) and for mid rapidity -0.5<y<0.5 

To check the consistency in our results we use two different partonic structure 
functions one of which is CTEQ5M and other an old one MRS125. The comparison 



-> — I — 1 — J—^ 

n F =n R =2Vm;+p T - 

m = 1.5 GeV 



1 1 1 1 1 1 

— a P t = (i bI , t> ci b = ,U2 

— Ap T = a L Vp T ,a L = 0.25 GeV 1 

— No Energy Loss 

p T < 2 GeV, lyl < 0.5 




Pb+Pb(8>LHC,Vs = 2.76ATeV 
_J i I i I i_ 



0.5 



1 



1.5 2 

A§ (rad) 



2.5 




3.5 



R 

ID- 



'S- m 3 

< io 

- -4 

Zio 

T3 



"i — 1 — i — 1 — i — 1 j ' ; i — 1 — i — 

, i 2 2 — No Energy Loss 

m =1.5 GeV - Ap T = 0,^,0^0.25 GeV 

_ Ap T = a B p T ,a B = 0.12 

p > 6 GeV, lyl < 0.5 




Pb+Pb@LHC,\'s = 2.76ATeV 
_i I i I i I i L 



05 



1.5 2 

Ad) (rad) 



2.5 



3.5 



Figure 3. (Colour on-line)same as Fig. 2, dN/dA(f> vs Acf> of DD pair for (left)py < 
2.0 GeV, (right)p T > 6.0 GeV. 
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Figure 4. (Colour on-line)Azimuhtal correlation of cc pair for (left) < pr >=2.0 GeV, 
(right) < p T >=3.0 GeV 



is shown in Fig. [H where the difference in the two distributions is very small and the 
shape almost identical. However more recent structure functions like CTEQ6M and 
CTEQ6.6 etc. must be used in order to have more up-to-date results. These issues will 
be addressed in our next publication on heavy quark correlation. 

Next let us recall that LO contribution can be differentiated from NLO contribution 
with different pt cuts on charm momentum. Leading order processes give back to 
back charm pairs which are entirely visible around A0=7r while NLO contribution is 
distributed from A0=O - it. 
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In Fig. El we show our results for dN c5 /d(j) for different px cuts. Realizing that 
all heavy quarks now appear with reduced momenta, we see that if we look at px < 2 
GeV or px < Q GeV, then the back-to-back correlation rise by up to a factor of 10 for 
= 0. The results for px > 2 GeV or px > 6 GeV are more dramatic in that the 0=7r 
correlation now reduces by more than a factor of 10 while that for 0=0 decreases from 
its value for no energy loss. 

We show dN D £)/dA(j) for p T > 6.0 GeV and px < 2.0 GeV in Fig. [3j Comparing 
it with Fig. [2] for same px regions, we observe certain differences which we now discuss. 
For px < 2.0 GeV, we observe that D meson distribution is slightly higher than charm 
spectrum at A0 = ir, although the order of magnitude remains same. While at A0 = 
0, the situation is reversed. Similar observations are noted when figures at px > 6.0 
GeV are compared. We feel that the above differences are caused by fragmentation 
function, D(z), which changes the px distribution of charm into px distribution of D 
mesons with, 0<z<l. Thus the correlation spectra of charm and D mesons may appear 
slightly different when we look into particular px regions. Finally it can be mentioned, 
D mesons rather than charms are observed in experiments. So calculating D meson 
correlation and comparing it with charm will give us deeper insights into the correlation 
study. 

In Fig. HJ we have E c E 5 dN/ d 3 p c d 3 p 5 for charm average px of 2 GeV and 3 GeV. 
The figure shows change in azimuthal correlation of charm pairs with pairs at A0=7r 
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decreased considerably by inclusion of the energy loss mechanism. 

To discuss our simple model of charm quark energy loss, we find that most of the 
charm pairs not only lose energy to shift to the lower momentum region but also back- 
to-back correlation for many charm pair is altered to almost collinear pairs. Also we 
find that two different energy loss mechanisms included in our study do not give much 
different outcomes. Further investigating at much higher momentum regions might bring 
out the differences between various energy loss mechanisms. The correlation study can 
be enriched if expanding medium is included in addition to energy loss by charms. 

Next we move to our results for correlated decay of charm. In Fig. [5j we have 
dN/dM e + e - for di-electrons from correlated charm decay. We can recall that there is 
enhancement in D mesons as well single non-photonic electrons due to the effects of 
large drag on charm quark moving through QGP. Here we find a similar enhancement 
in di-electron spectrum at midrapidity. For M e + e - less than 1 GeV, there is increase 
in dN / dM e + e - by almost 12% which is quite noteworthy considering our model to be 
simple empirical mechanism of energy loss. 

5. Summary 

We have studied correlation of charm, D mesons as well as correlated decay of charm 
using NLO-pQCD processes and a simple empirical model for energy loss. The azimuthal 
correlations of charm show change when energy loss mechanisms are implemented along 
with cuts on charm transverse momentum. In case of di-electron distribution, energy 
loss enhances the electron spectrum slightly for low invariant mass. 
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